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Abstract. The observed mass-to-light (M/L) ratios of globular clusters 
(GCs) are on average ~ 20% lower than expected from Simple Stellar Population 
(SSP) models, which only account for the effects of stellar evolution. We study 
the M/L ratio evolution of a sample of 24 Galactic GCs using parameterised 
cluster models. The dynamical evolution of GCs is included by accounting for 
their dissolution and by using a detailed description of the evolution of the stellar 
mass function. The ejection of low- mass stars leads to a decrease of M/L, which 
is found to explain the discrepancy between the observations and SSP models. 



1 Introduction 

Galaxy mergers are thought to have been the formation sites of globular clus- 
ter systems (GCSs). As such, globular clusters (GCs) can be used as trac- 
ers for studying the early evolution of galaxies and any mergers these galaxies 
may have experienced. To relate the GCS of a galaxy to its formation history, 
it is essential to obtain a good consensus on the present masses of the GCs. 
For obtaining these , (constant) mass-to-light (M/L) ratios are c ommonly used 
(|Fall fc Zhandl200ll : [Jordan et alJl2007l : lMcLaughlin fc FaJil200i ). However, the 
observed dynamical M/L ratios of GCs in several galaxies are found to be ~ 20% 
low er than the values exp e cted from Simple Stellar Population (SSP) mod- 
els (iMandushev et al.lll99ll ; iMcLaughlin van der Marell [20051 ; IReikuba etH 
20071 ). This complicates the interpretation of GC masses, as it first needs to be 
understood why their M/L ratios deviate from those predicted by SSP models^. 

It has been pr oposed that the M/L ratio difference is du e to the dynamical 
evolution of GCs (|Kruijssenll2008l ; Ikruiissen fc Lameri 120081 ) . The preferential 
ejection of low-mass, high-M/L stars from dissolving star clusters changes the 
shape of the stellar mass function (MF) within a cluster ( Vesperini fc Heggid 



Interestingly, ultra-compact dwarf galaxies (UCDs) , which have been proposed to represent 
a continuation of the GC mass range t o higher masses, have M/L ratios that ar e ~ 25% higher 
than those predicted by SSP models (IReikuba et al.1 120071 ; iMieske et "all 120081 ). We separate 
this discrepancy from the one between GCs and SSP models and only consider GCs. 
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19971 : iBaumgardt k Makinol [20031 ; |Pe Marchi et al.ll2007l ). causing the luminos- 
ity to be only marginally affected while the mass decreases. As such, dissolution 
decreases the M/L ratio with respect to the expected evolution from SSP mod- 
els, in which it is assumed t hat th e shape of the stellar MF does not vary. 

In [Kruiissen &: Mieskd (|2009l ) , the hypothesis of low-mass star depletion as 
an explanation for the low M/L ratios of GCs was tested by consid ering the sam- 



ple of 24 Galactic GCs for which the orbits (Dinescu et al1ll999l ) and observed 
dynamical M/L ratios ( McLaughlin &: van der Marelll2005l ) are known. There, 



we derived the dissolution timescales from the individual cluster orbits and com- 
puted the resulting M/L evolution for the GCs in the sample. It was found that 
low-mass star depletion can indeed account for the ~ 20% gap between the 
observations and SSP models. As a first-order approximation, low-mass star 
depletion was included by increasing the lower stellar mass limit in the cluster, 
which was tuned to iV-body simulations in order to give reliable results. In 
reality, the stellar MF evolves more gradually. This has been included in new, 
physical models of the evolution of the stellar MF in dissolving star clusters 
( Kruiissen 2009). The aim of the present paper is to revisit the calculations of 
Kruiissen Sz Mieskd ( 20091 ) with these new models and to verify whether their 
conclusions still hold. 



2 Star cluster evolution and mass-to-light ratio 



We u se the parameterised cluster model SPACE (IKruiissen fc Lamersl[2008l : lKruii sscn 
2009), which incorporates the effects of stellar evolution, stellar remnant produc- 
tion, dynamical dissol ution and energy equipartiti on. The dissolution timescale 



due t o evaporation (cf. IBaumgardt Sz Makinoll2003l ) and tidal shocks (cf. lDinescu et al 
1999) is determin ed for each of the 24 GCs in our sample by considering their 
individual orbi ts (IKruiissen Sz Mieske 20091 ) and is subsequently converted to a 
mass loss rate (jLamers et al.ll2005l ). 

The evolution of the stel lar MF is comp uted by considering the ejection rate 
as a function of stellar mass ( Kruiissen 20091 ). The adopted method accounts for 
mass segregation and dissolution in a tidal field by using the timescale on which 
energy equipartition is reached for different stellar masses and by comparing 
the stellar velocities with the escape velocity. The p hotometry is comp uted by 
integrating the MF over the new Padova isochrones ( Marigo et al.1120081 ). 



3 Comparison to observations 

The y-band M/L ratios from SSP models are compared to the observed values 
in the left-hand panel of Fig. [TJ In the right-hand panel, the modeled y-band 
M/L ratios for the 24 GCs in our sample are compared to the observations. 
The improvement with respect to the left-hand panel is considerable, but there 
remains a number of GCs with obse rved M/L ratios th a t are lower than the 
modeled ones. This was explained by IKruiissen Sz Mieske ( 2009), who reasoned 
that some of the observed M/L ratios from lMcLaughlin van der Marell (|2005l 1 
are likely biased due to their use of isotropic single-mass King models to deter- 
mine the M/L ratios. This method leads to a tendency towards the central M/L 
ratio rather than the global one for GCs with evolved internal structure. 
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Figure 1. Comparison of F-band M/L ratios that are predicted by SSP 
models (left) with observed M/Ly ratios, and of our modeled M/Ly ratios 
(right) with the observed values. Error bars denote lcr standard errors, with 
diamonds marking GCs for which theory and observation are in la agreement, 
and dots representing the remaining deviant GCs. 



In total, our new models explain the M/L ratios of 17 out of 24 GCs 
wi thin the lcr error margins . This is a substantial development with respect 
to Kruiissen Sz Mieske (2009), where half of the GCs was explained. The use of 
the new MF models and of the new Padova isochrones contribute equally to this 
improvement (three GCs each). The average fraction of the observed M/Ly ra- 
tio with respect to the value predicted by SSP models is 0.81^QQg, while for our 
modeled fraction this is 0.79±0.01. These values are in excellent agreement. An 
independent check of our mode ls is provided by compa ring the modeled stellar 
MF slopes to those observed bv lDe Marchi et all ( 20071 ). This is shown in Fig. [2] 
and confirms the general validity of our models, despite the scatter and large 
error bars. 



4 Discussion 



Low-mass star depletion in dissolving GCs explains the ~ 20% discrepancy be- 
tween the obse rved M/L ratios and those predicted by SSP models. The new 
cluster models ( Kruiisse"nll2009l ). in which we account for the changing slope of 
the stellar MF rather than shift i ng th e lower stellar mass limit, show that the 
results from lKruiissen & Mieskel (|2009l ) also hold when more detailed models are 
applied, and improve the agreement between theory and observations. 

Nonetheless, care should be taken when using observed dynamical M/L 
ratios that neglect a stellar mass spectrum or the variability of M/L from the 
centre to the outskirts of a cluster. These cannot be interpreted as global M/L 
ratios if the M/L too strongly varies of radius. Considering the central role GCs 
play in (extra) galactic astronomy, it essential to obtain an accurate census of 
their dynamical masses. This would improve the observed globular cluster mass 
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Figure 2. Observed stellar MF slope a bs for the stellar mass range 0.3 < 
to/Mq < 0.8 versus the modeled mean slope in that range a pr od. Error bars 
are la standard errors, mainly caused by the uncertainty in the GC mass. 



function ( Kruijssen Sz Portegies Zwartll2009l ) . which can then be more accurately 



used to trace and interpret the formation history of galaxies. 
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